400 



Analysis of Packing Function 
Solutions for Monomelic 
Proteins 




9 12 15 18 21 24 zTso 
Distance (A) from 
Experimental Average 
Coordinate 



500 




% 4 



Distance {A) from 
Experimental Average 
Coordinate 



24 27 



START: 



SctUp 
Parallelization 



Initialize 
Constants and 
Sine/Cosine Tables 



For Use with a Radial Correction or 
with Modes 5 thru 7: 



Input 1 & 
Initialize the List of 
Fourier Indices (hkl) 



Get Biggest Sphere 
Without Overlap 
of Symmetry Mates 



Get 

Fractionaiizalion 
Matrix 



Partition the List 
of Fourier Indices 
Between 
Proccssorfs) 



Convert Each Fourier 
Index (hkl) to Polar 
Coordinates & Get 

IU Bcsscl Argument 



Confirm 
Calculation Mode 



Get the FT of 
Crystalline Unit Ceil 
Filled with Symmetry 
Related Spheres 



Calculate a Radial 
Correction Factor 
for Each "n" Index 



Initialize 
Tables of Spherical 
Bessel Functions and 
Besscl Function O's 



Determine Limits 
on 1 & m indices and 
on n indices at each 1 



For modes 5 thru 7: 
Read expected SHSB 
coefficient & 
deviation values 



Calculation Mode-Specific Routines: 

Modes I & 2 
(Unphased Diffraction Amplitudes to 
Phased FT ofSHSB-modeled Unit Cell) 



Mode-Specific 
Input 



Get Input File # for 
Calculated FT, if 
there is a Prior Model 



Initialize Correlation 
Exponentiation 
Factor 



?low Chan for die Main Drivcx Program for "fiuzer": 



Modes 1& 2 (cont'd) 



Cycle 

For Each Choice (mstop) of 
Stopping Value for the SHSB 
m index 



Update Correlation 
Exponentiation 
Factor 



If there Is a Prior Model ™ 
(On cycles after the 1st cycle 
of the 1st run of the program): 



On Cycles After the 
1 st value of mstop 
Update the File # for 
the Calculated FT 



For Each input SHSB Model 



Convert the File # 
& the Model # 
to a File Name 



Read Input 
FT & SHSB coefs 
of the Prior Model 



If" here Is No Prior Model: T 



Model the Crystal as 
a Crystalline Unit Cell 
Filled with Symmetry 
Related Spheres 



Get a Starting Value I 
(msus) for the SHSB I 
m Index for the Next! 
Packet of m Values I 



Update the Model 
by Adjusting the 
Contributions of 
each SHSB function 
for each lmn index in 
the m-packet from 
'msus' to 'mstop' 



': Optiona to compote a the FT of * SHSB Model 



of Crystal 



Initialize 
Fractlonalization 
Matrix 



initialize the Equal 
Partitioning of the 
Fourier (hkl) Index 
between Processors 



On 1 st Cycle of 1 st Run ; 



Prescale Observed 
Diffraction to that 
of a Unit Cell of 
Spheres 



Define the First 
SHSB Index Triplet 
(Imn) for which to 
Consider Model |F|'s 



Initialize for Index- 
by-lndex Update 
of Origin-Centered 
SHSB Basis Function 



Modes 4 & 5 only: 
Initialize Buffers for 
Cumulative Update of 
Fourier Representation 



Initialize Pointers to 

Stored Fourier 
Representations of 
Model and of Basis 



Mode 3 only; 
Get File Name from 
File #& Open It to Let 
5H5BCoefs.be Read 



For each "m" Index 
(0 to maximum *m" ) 



For Each hkl in this 

Partition: 
Update "m" Recursion 

Formula for Fourier 
Representation of the 
Origin-Centered SHSB 



For each "I" Index 
(present "m* to maximum T) 



For Each hkl in this 
Partition: 

Update "l,m" Recursion 

Formula for Fourier 
Representation of the 
Origin-Centered SHSB 



For each "n" Index 
(1 to maximum "n" for each "I*) 



For Each hkl in this 

Partition: 
Update "n" Recursion 

Formula for Fourier 
Representation of the 
Origin-Centered SHSB 



Depending on Mode: 
Choose the * of Passes 

and # of Presumed 
Phase Angles Needed 
for the SHSB coef. with 
this SHSB index (Imn) 



Set the Presumed 
Amplitude of the 
Origin-Centered 
SHSB Basis Function 



FIRST PASS: 

initialize Registers: 
Overall Comparison 
of Correlation Coef. 
& Other Statistics 



Renitialize Pointers to 
Storage Sites for Fourier 

Representations of the 
Full-Unit-Cell SHSB Basis 



Parallel Processor Version: 



Set # of calculations to: 
(# of presumed values of 
SHSB coef. "s phase) 
X 

(# of stored accumulated 
SHSB models for trial 
combination with this 
new SHSB component) 



Given: # of processors 

# of hkl partitions 

# of calculations 
Get* of required 

rounds of trial 
combinations 



I For each round of trial 
combination on this processor 



Single Processor Version: 
(Outer Loop) 



For each presumed value 
of the SHSB coef.'s phase 



initialize Registers: 
Angular Comparison 
of Correlation Coef. 
& Other Statistics 



